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Abstract 



We calculate the trade-off between the quality of estimating the quantum 
state of an ensemble of identically prepared qubits and the minimum level 
of disturbance that has to be introduced by this procedure in quantum me- 
chanics. The trade-off is quantified using two mean fidelities: the operation 
fidelity which characterizes the average resemblance of the final qubit state 
to the initial one, and the estimation fidelity describing the quality of the ob- 
tained estimate. We analyze properties of quantum operations saturating the 
achievability bound for the operation fidelity versus the estimation fidelity, 
which allows us to reduce substantially the complexity of the problem of find- 
ing the trade-off curve. The reduced optimization problem has the form of 
an eigenvalue problem for a set of tridiagonal matrices, and it can be easily 
solved using standard numerical tools. 
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I. INTRODUCTION 



It is a well-known fact that given a single copy of a quantum system it is in general 
impossible to determine its quantum state exactly This principle is closely related to 
the no-cloning theorem || , which prevents one from producing multiple faithful copies from 
a single unknown state. Of course, the situation becomes different when one is given an 
ensemble of identically prepared quantum systems. With the increasing size of such an 
ensemble one can extract more and more precise information on its preparation M. 

An important effect which usually accompanies an attempt to determine the unknown 
preparation of a quantum system is the disturbance of its original state. The state distur- 
bance is a penalty for gaining classical information from a quantum system. This fundamen- 
tal feature of quantum mechanics has been discussed from many different points of view, 
depending on the particular physical scenario considered |4|,[| . A recent paper || presented 
the description of the trade-off between the information gain and the quantum state distur- 
bance motivated by the problem of quantum state estimation [||,[7|-[I(J . In this approach, the 
classical outcome gained from an operation on the quantum system is converted into a guess 
what the original state was. According to the general quantum mechanical rule linking the 
information gain with the state disturbance, the better guess can be made on average, the 
less the final state of the system should resemble the initial one. Natural parameters for 
quantifying the trade-off in this context are mean fidelities, defined using the scalar product 
between the relevant state vectors, and averaged over many realizations of the scenario. 
The first quantity of interest is the operation fidelity F, which parameterizes the average 
resemblance of the state of the system after the operation to the original one. The second 
quantity, known as the estimation fidelity G, tells us how good estimate one can provide on 
the basis of the classical outcome of a given operation. In the plane of the fidelities F and G, 
quantum mechanics imposes a bound in the form of a trade-off curve limiting the maximum 
values of F and G that can be achieved simultaneously for any quantum operation. For a 
single copy of a d-level system, the trade-off curve between F and G can be described in 
simple analytical terms 0. 

In this paper we discuss a more general case of the fidelity trade-off, when one is given 
a finite ensemble of identically prepared systems. We shall assume that these systems are 
qubits prepared in an arbitrary, randomly selected, pure state. Of course, the penalty — i.e. 
the state disturbance — for gaining information is expected to be smaller for an ensemble 
compared to the single system case. In the limit of an infinite number of copies, we can 
determine the quantum state exactly and reset the state of all the qubits according to this 
precisely known information. Then in principle no disturbance needs to occur. Our interest 
here will be focused on the intermediate case between the single-copy operations and nearly 
perfect estimation of large ensembles. We shall discuss the trade-off curve between the 
operation fidelity and the estimation fidelity in the most general case when the ensemble 
consists of a finite number N of qubits. 

The two fidelities used in our work have fundamentally different practical meanings. 
Informally speaking, the operation fidelity F deals with the intrinsic quantum information 
remaining in the state of a quantum system. In contrast, the estimation fidelity G describes 
the classical information gained from the measurement. This classical information allows us, 
for example, to generate arbitrarily many copies of the initial state with the same fidelity G. 
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Consequently, for finite ensembles the fidelities F and G take values from different ranges. 
For example F can easily be equal to one for any size of the ensemble (which simply means 
that nothing is done to the qubits) but in contrast G = 1 means that we are able to generate 
arbitrarily many perfect copies using the result of the estimation. Previous work on quantum 
state estimation H demonstrated that the maximum attainable estimation fidelity for an 
ensemble consisting of iV qubits is G = (N+1)/(N + 2). As we justify later, this is also the 
value for the operation fidelity in the limit of optimal quantum state estimation. Of interest 
here is what happens below this value for G and above this value for F, as this describes 
the region where we try to convert part of the the quantum information contained in the 
initial state into a classical guess. 

Finding the trade-off curve for an ensemble of N qubits presents in principle a rather 
complicated problem. A general quantum operation can map the initial state of the qubits 
onto the full Hilbert space with the dimensionality 2^, which grows exponentially with the 
size of the ensemble. Furthermore, the classical operation outcomes used for the estimation 
can a priori assume values from an arbitrarily large set. In this paper we demonstrate in 
several steps that it is possible to reduce the general problem of finding the trade-off curve to 
0(y~N) independent constrained optimization problems, each involving O(N) real variables. 
The reduced optimization problems have quadratic form, and they can be solved numerically 
by finding the eigenvectors of certain tridiagonal matrices. This is a substantial reduction 
of the complexity of the problem compared to its original formulation, which allows us to 
deal numerically with much larger ensembles. Furthermore, numerical evidence strongly 
suggests that in general just one from 0{yN) optimization problems gives the full trade-off 
curve, but a strict mathematical proof of this conjecture is lacking. Most importantly, we 
achieve the reduction of the complexity without imposing any restrictions on the generality 
of quantum operations considered, and the calculated trade-off curves are both tight and 
universal. 

The results of our paper are summarized in Fig. |Tj, where we depict the trade-off curves 
calculated using our approach for several exemplary values of N . All the curves have a 
common extreme point attained for F = 1 and G — 1/2. This point is reached when the 
ensemble is simply left intact, and the corresponding value for the estimation fidelity G = 1/2 
describes making a completely random guess about the state of the qubits. The other 
extreme point for each trade-off curve, corresponding to optimal quantum state estimation, 
is given by F = G = (JV + 1)/(N + 2). In this limit, the fidelity of the qubits remaining 
after the operation is exactly the same as the fidelity of our guess. One can give a simple 
intuitive argument that this should be the case: F is always an upper bound on G, since 
one can always set the state of the qubits equal to the guess state. It is plausible that 
this bound can be achieved in the limit where we only care about maximizing G. Between 
the two extreme points, the depicted trade-off curves illustrate how with increasing N the 
extraction of classical information has less of an effect on the state of the qubits after the 
operation. 

The problem considered in this paper can be viewed as a special case of quantum cloning, 



i.e. generating a larger number of imperfect copies from a given ensemble [11]. As we noted, 



given a classical estimate of the quantum state, we can use it to generate an infinite number 



of qubits with the same fidelity as the estimate [12|. Thus, the trade-off curves presented 



in this paper describe the optimal performance of an asymmetric quantum copying machine 



3 



which given N identical pure qubits produces N clones with the fidelity F, and in addition 
arbitrarily many clones with the fidelity G. 

The paper is organized as follows. First, in Sec. [TIL we formulate the problem of the 
fidelity trade-off in quantitative terms. In Sec. |T| we simplify the formulae for the fidelities 
using the angular momentum representation of the rotation group. This provides explicit 
expressions for the fidelities that are suitable for further calculations. In Sec. [IV| we argue 
that in order to find the trade-off curve it suffices to consider a single quantum operation 
element, thus substantially reducing the complexity of the problem. We further demonstrate 
in Sec. [V] that it is sufficient to consider operations which map the initial state of the 
qubits only onto the fully symmetric subspace. With these results in hand, we define in 
Sec. |VI] the reduced optimization problem which yields the actual trade-off curve. The 
numerical solution to this problem is discussed in Sec. |V11| . Next, we show in Sec. |V111] that 
the calculated trade-off curves are achievable, by constructing explicit quantum operations 
attaining the derived bound. Sec. [TX| concludes the paper. 



II. FORMULATION OF THE PROBLEM 

We begin with an ensemble of N qubits all prepared in the same pure state |0). We 
shall use the following notation: 

|n) = #(n)|t>, (2.i) 

i.e. the state |0) is represented as a result of a unitary operation U(Q) on a reference state 
| j), which we will take for concreteness to be the spin up state along the z axis, a z \ j) = | f )• 
The group of the unitary transformations U(Q) can be conveniently parameterized using the 
two-dimensional irreducible representation of the rotation group. Thus Q can be considered 
as an abbreviation for the triplet of the Euler angles (</>, 6, (), and 

U(Q) = exp(-i(j)a z /2) exp(-t6a y /2) exp(-i(a z /2). (2.2) 

The third Euler angle ( introduces a trivial overall phase factor in the definition of the 
states \ and in principle it could be set to zero. However we will keep it as an independent 
variable in order to adhere strictly to standard angular momentum algebra notation which 
we will use later. The canonical volume element in the group of unitary transformations 
U(fl) is given by 



dQ = ^sm9d9d(pdC. (2.3) 



We assume that the initial state of the ensemble of N qubits is randomly selected 
according to the probability distribution given by this measure. 

Initially, the composite state of the ensemble of the qubits is described by a tensor product 
|r2)(fi|® Ar . We assume that the qubits are submitted to an action of a certain quantum 
operation, which can in general act collectively on the whole ensemble. Such an operation 
is described by a set of operators {v4 rs } acting in the 2 7V -dimensional tensor product Hilbert 



space of all the qubits |0. The classical outcome of the operation is given by the index 
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r, and it is correlated with the final quantum state of the qubits. The probability p r (0) of 
obtaining the result r is given by: 

Pr (n) = J2M^rsArM(nr N ). (2.4) 

s 

The conditional transformation of the ensemble corresponding to the outcome r is described 
by the formula: 

g*(n) = -^£i^)^riL- (2-5) 

The summation over the index s maps in general pure states onto mixed ones, and it can be 



viewed as responsible for introducing excess stochastic fluctuations |TJ[]. In general, the index 
s can assume values from a different set for each r. For the operation to be trace-preserving, 
the set of the operators {A rs } must satisfy the completeness relation of the form: 

£it s i rg = l (2 .6) 

rs 

With the above notation for quantum operations, we can now define explicitly the two 
quantities central to this paper: the operation fidelity F and the estimation fidelity G. The 
operation fidelity F quantifies the average resemblance of the state after the operation to 
the original one. Let us consider the reduced single-qubit density matrix after the operation, 
averaged over all N qubits: 

Ql cd ({l) = iflW^n) + Tr 1)3 -iv£ ut (^) + • • • + Tr HN -i)Q° r nt mi (2-7) 

where the subscript of the trace symbol Tr labels the range of qubits over which the trace 
operation is performed. The expectation value of the above expression on the initial state 
(fi|^ cd (fi)|fi) tells us how much on average the state of a single qubit after the operation 
resembles its initial value. This quantity, summed over the possible outcomes r of the 
operation with the corresponding weights p r and averaged over the initial state of the qubits, 
yields the mean operation fidelity for an ensemble of identically prepared qubits: 

F= /dfi£p r (ft)(fi|g cd (ft)|Q). (2.8) 
J r 

The second quantity of interest is the estimation fidelity G. Given the classical outcome 
r of the operation, we can make a guess \Q r ) what the original state of the qubits was. A 
natural way to quantify the quality of the guess is to take the squared absolute value of the 
scalar product between the guess and the original state, equal to |(fi r |fi)| 2 . The estimation 
fidelity is obtained by averaging this expression over the sets of possible operation outcomes 
r and the input states \ fl): 

G= f dQ^2p r (Q)\(Q r \Q)\ 2 . (2.9) 

r 

The estimation fidelity depends not only on the quantum operation {A rs } itself, but also 
on the estimation rule used to make the guess, described by the mapping r i— > \Q r ). We 
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will demonstrate in the following how to define this mapping in a way which optimizes the 
estimation fidelity for a given arbitrary quantum operation. 

Our goal is now to find the inequality which bounds the fidelities F and G, assuming the 
most general form of the operation {A rs }. The only requirement that has to be satisfied by 
the set of the operators A rs is the completeness condition described by Eq. (|2.6| ). 



III. FIDELITIES 

In this section, we will simplify the expression for the fidelities F and G to the form which 
makes them more manageable in the optimization procedure. Our first step will be explicit 
calculation of the integrals over the space of pure states \Q), which can be done with the help 
of tools developed in the theory of representations of the rotation group. Throughout this 
section, we shall follow strictly the notation of Ref. |y| for the angular momentum algebra 
and the elements of rotation matrices. 

To simplify subsequent expressions, we begin with a general observation that the index 
s appearing in Eqs. ( |2.4j ) and ( |2.5| ) could in principle also be known classically after the 
operation. Summation over the index s in Eqs. (|2.4| ) and Eq. ( |2.5| ) means that the operation 



is imperfect, and it averages statistically different output states. This results in the loss of a 
fraction of the information extracted from the initial quantum state. As we are interested in 
the optimal operations saturating the quantum-mechanical bound on the fidelities, we can 
assume with no loss of generality that both r and s are known. In such a case, we can use a 
new single index to label both r and s. For this reason, we will assume in the following that 
the index s is trivial, i.e. it assumes only a single value for each r, and that consequently it 
can be dropped from further notation. Thus we restrict our attention to quantum operations 
which are known in the literature as ideal [[16] or efficient |J. 

In the following calculations, it will be convenient to use the decomposition of the com- 
plete Hilbert space of N qubits into subspaces with the fixed value of the total angular 



momentum operator. This decomposition has the form |T7 



®n=®n Ja , (3.i) 

i=l a=l 

where j = N/2 is the largest value of the angular momentum appearing in the decomposition, 
the j^s assume values in the set {j,j — 1, ... ,j — \J\} and are arranged in a non-ascending 
order. Here 7ij> denotes the subspace corresponding to the total angular momentum value 
j'. We note that j\ — 3 an d j a < j for a > 1, as there is only one representation with 
j = N/2 corresponding to the fully symmetric subspace. For completeness, in Appendix |A| 
we present a simple derivation of the multiplicities of the angular momentum representations 
appearing in Eq. (|3.1|) . The action of a tensor product of the unitary operators U® N (Q) in 
each of the subspaces Hf is given by the corresponding representation of the rotation group. 
In order to simplify the notation, we will henceforth use the same symbol U(Q) to denote 
the action of U® N (Q) on the whole ensembles of qubits, as for a single qubit. 
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A. Operation fidelity 

Let us start with the operation fidelity F. The matrix element (Q\gl° d (Q)\Q) appearing 
in the definition of F in Eq. ( |2.8| ) is equivalently given by the expectation value over the 
density matrix g° nt (Q) defined in Eq. (|2.5| ) of the following operator: 

P{Q) = + \Q)(Q\ ®...®i + i(8)i(8)...(8) |0)(0|) (3.2) 

Using the above definition, the mean estimation fidelity can be written compactly as: 

F = [ dn^Tr[P(p)Ar\n)(n\ m Al]. (3.3) 

We will now express P(fl) in terms of the angular momentum operators. Expressing the 
projection |f2)(f2| as a combination of the Pauli matrices: 

\n)(n\ = ^(i + sin9cos0a x + sin9sin(pa y + cos9a z ) (3.4) 

we can write the operator P(Q) in terms of the total angular momentum operators for the 
composite system of N qubits: 

1 N . 

P(Q) = — - ^(1 + sin 9 cos 4> of + sin 9 sin a\ + cos 9 of) 
^ i=i 

= + ^(sin6»cos0J x ' + sin 9 sin 4>J y + cos9J z ), (3.5) 

where the index i enumerates the qubits. In the following, it will be convenient to switch 
to the pair of the angular momentum raising and lowering operators = J x ±iJ y , using 
which we have: 

P(O) = -U + — (~e*+ sin 9 J- + -e* sin 9 J + + cos 6 J*} (3.6) 



With this expression for the operator P(fi), we can write the mean fidelity F as: 
F = \ + ^ E Tr (^J-AK^Al + -j+A^Al + J z A r K A^j , 



(3.7) 



where all the terms involving the Euler angles Q have been collected to three integrals K T , 
r = —1, 0, 1, defined as 

K T = J dttk T {tt)\tt)(tt\® N (3.8) 

with the functions k T (Q) given by: 

fc_i(fi) = e^sinfl 
k (Q) = cos 9 

k^Q) = e-^smO. (3.9) 
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Explicit calculation of the integrals K T , performed in Appendix [B[ yields the following 
expressions: 

1 A=F 
K±1 = (j + 1)(2j + 1)^ 

^-(TTiyWTi)^ (3 - 10) 

where j = iV/2 and 3f , Jj denote the angular momentum operators restricted to the com- 
pletely symmetric subspace of the N qubits defined by the angular momentum j = N/2. The 
operators K T vanish outside this space, as they are defined as integrals of totally symmetric 
projection operators \ fl) (fl\® N . 

Inserting the explicit form of the operators K T into Eq. ( ft.7|) yields: 

F = \ + w + m + 1) ? Tr (¥~ AJ ? Al + l rU i! + J ' AJ ' Al ) (3 - n) 

The operation fidelity can be equivalently expressed in terms of the Hermitian operators J x , 
J y , and J z as: 

F = \ + 2 j(j + l)(2j + l) ? ' VnrA ^^ + J VA ^ A 1 + J Z AJ!At) (3.12) 

We note that this expression is completely symmetric with respect to the three Cartesian 
components of the angular momentum. Furthermore, it can be easily checked that each of 
the traces appearing in the sum over r in the above formula is invariant with respect to the 
transformation of the operator A r according to: 

A r -> u(n)A r fr(n), (3.13) 

where U(Q) is an arbitrary rotation matrix. 



B. Estimation fidelity 

We will now evaluate the integral over the space of pure states in the expression for the 
estimation fidelity G, given by: 

G= [ dnJ2^( A iAr\ty(n\® N )\(n r \tt}\ 2 (3.14) 

J r 

As we discuss in Appendix |B|, the projection on the product state |f2)(f2| (g,Ar can be repre- 
sented as: 

\n)(n\® N = u(n)\j; 3 )(j;j\u\n), (3.15) 

where j = N/2 and \j;j) belongs to the subspace with the total angular momentum j and 
is the eigenvector of J z corresponding to the eigenvalue j. Thus we have: 
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(3.16) 



We can now swap the order of the integration over Q and the summation over r, and change 
for each r the integration variables from Q to Q' such that: 

u(n') = u\n r )u(n) (3.17) 

Using this parameterization, we have (f |?7^(f2 r )f/(f2)| |)| 2 = (1 + cos#')/2 and 

G = \^j ^'(I + cos^^^W^OI^^O'JJI^W^W] (3.18) 

This expression can be decomposed into two parts according to the two terms in the factor 
1 + cos^'. The first part involves the integral 



d^U(^)\j;j)(j;j\U\n') 



2J + 1 J 



(3.19) 



which is proportional to the identity operator truncated to the completely symmetric 
subspace. Consequently, the summation over r for this term can be easily performed which 
yields the constant value 1/2. The second part can be expressed with the help of the operator 
Kq defined in Eq. (|3.8| ), which gives: 



(3.20) 



Q = 1 ^ Tr(AtA r ^) + 1 Tr[i ti^ (a) ^ o ^t (a)] 



Inserting the explicit form of K yields: 



i 



2(j + l)(2j + l) 



^Tr[iti r ?/(fi r )jjf>t( fir )] 



(3.21) 



This expression for the estimation fidelity allows one to derive easily the optimal estimation 
strategy, i.e. the mapping from the set of outcomes r to guesses Q r for a given quantum 
operation {A r }. In order to derive this strategy we first note that 



U(Q r )JjW(n r ) = sin 9 r cos 4> r JJ + sin 9 r sin <p r Jj + cos 9 r J?. 



(3.22) 



This allows us to write each of the traces in the sum over r in the form of a scalar product 
AiTfL between two three-dimensional real vectors A r and Q r defined as: 



/ Tr(iti r J?) \ 

Tr (it i r JJ) 
\ Tr(iti r j|) J 



(3.23) 



and 




(3.24) 
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It is seen that Q r is a unit vector pointing in the direction defined by the Euler angles Q r . 
Obviously, the scalar product AjQ r is maximized when the two vectors are parallel. Thus, 
for a given r the corresponding Euler angles <p r and 9 r should be chosen such that: 




(3.25) 



Of course, the value of the third Euler angle ( r can be arbitrary, as it introduces only an 
irrelevant overall phase factor. If the vector A r is zero, then the guess has to be made 
completely randomly. Thus, the maximum value of the trace for a specific r in Eq. (|3.21| ) is 
given by \J AjA r , and the maximum value of the estimation fidelity attainable for a given 
operation {A r } is defined by: 

G = \+ 2{J + 1) 1 (2j + 1} E TjmttArJ?)]* + {Tr(MA r J])}i + [Tr(MA r J*W. (3.26) 

As in the case of the operation fidelity, this expression is invariant with respect to an arbitrary 
transformation of the operators A r according to the rotation group. 



IV. DECOMPOSITION 

Let us now summarize the constrained optimization problem describing the trade-off 
between the operation fidelity F and the estimation fidelity G. These two quantities can be 
written as: 

F = - + —^ -f (4.1) 



and 



G = - + , 1 , g (4.2) 
2 2(j + l) y 1 ; 



where the terms / and g depending explicitly on the operation {A r } are given by: 

1 

2j + 

and 



/ = 7-^- V Tr(J x A r J*Al + J y A r jyAl + J z A r .hA\) (4.3) 
2? + 1 ~f j j j 



9 = 2~~Y E \f[H^W + [Tr(iU r JJ)] 2 + [Tr(MA r J*f (4.4) 

It is worthwhile to look first at the extreme cases. The maximum value of the operation 
fidelity itself is of course F — 1. This corresponds to / = j(j + 1). This limit is achieved 
by the identity operation, for which it is easy to check that g = 0. Hence we can assume in 
the following that g is a nonnegative quantity. The other extreme case is the optimization 
of the estimation fidelity alone, which has been studied previously [0]. According to these 
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results, in the limit of optimal quantum estimation we obtain g = j, which sets the upper 
bound on the region of interest for g. We will demonstrate in Sec. |VII| that the maximum 
value of / attainable in the case of optimal estimation is equal to / = j 2 . Expressing this 
in terms of the fidelities, we get that F = G. 

The operators A r must satisfy in general the completeness condition described in 
Eq. (|2.6|) . However, since the initial state of the N qubits lies in the completely symmetric 
subspace, we need to consider the action of these operators only on the fully symmetric 
subspace described by Tij. Consequently for the purpose of our discussion we shall assume 
that A r : TCj — > TC® , and that the trace preserving condition takes the form: 

Y,%K = U (4-5) 

r 

Henceforth we replace the trace preserving condition by the milder requirement, obtained 
by taking the trace of the above equation: 

Tr^i^ = 2j + l. (4.6) 

We will perform the optimization under this weakened constraint and then show that the 
optimal (/, g) curve may be attained by a set of A r that is actually trace preserving. 

Weakening the completeness condition allows us to introduce an important simplification 
in further calculations. As we noted in Sec. |III A| , the expression for the operation fidelity 
is invariant with respect to rotations performed on the operators A r . This is also the 
case of the weaker trace preserving condition described in Eq. (|4.6j ). Consequently, we can 
always modify each of the operators A r by an operation of the form A r — > U(Q)A r W(Q) 
such that the vector defined in Eq. (|3.23| ) is aligned along the z axis and moreover its z 
component is nonnegative. In explicit terms, we assume that Tr(AJAjJ) = Tr(AJAjJ) = 0, 
and Tt(AIAJj) > 0. This allows us to replace the square root appearing in Eq. (|4.4|) by a 
much simpler expression Tr(AjA r JJ). An important advantage of this step is that the latter 

expression is quadratic in the matrix elements of the operators A r . 

Our next step will be the representation of / and g as linear combinations: 

f = J2Kf(A r ) 

r 

g = Y / Kg(A r ) (4.7) 

r 

where 

f(A) = L^Ti(J x AJ\ c i t + J y AJ y A ] + j z AJ*A r ) (4.8a) 

Tr(AU) 3 3 3 

g(A) = \ — — Tr(A^AJf) (4.8b) 

and the nonnegative coefficients A r are given by 

Trf^i) 

2j + l • { ' 
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The values of (f(A r ), g(A r )) are insensitive to the rescaling of A r , so A r are free variables 
up to the constraint 

E A r = l (4.10) 

r 

resulting from the weaker trace preserving condition. Hence each (/, g) point is a convex 
combination of independent (f(A r ), g(A r )). As we are interested in the boundary of allowed 
(/, g) points, it suffices to examine the case when the A r are all zero except for one value 
of r. This means that our problem is solved by using only one operation element which we 
denote by A, obeying the constraint Tr(iU) = j + 1. This observation has also been used 
in quantum rate-distortion theory ||18|| . After deriving the bound for (/, g) based on a single 
element A, we will demonstrate in Sec. [VII1| that the operator A can be used in a canonical 
way to construct a quantum operation satisfying the original full completeness condition. 



V. FULLY SYMMETRIC SUBSPACE 

We will now show that the optimization problem can be simplified even further: namely, 
that it is sufficient to consider operators A which do not transfer the state of N qubits 
beyond the fully symmetric subspace. 

According to our discussion of the structure of the Hilbert space of the whole ensemble, 
the operator A : Tlj — > TL® N may be viewed as consisting of entries A a , each entry acting 
from A a : Tlj — > TCj a - As the angular momentum operators are block diagonal and they do 
not mix subspaces with different a's, we can introduce the following decomposition of the 
quantities / and g: 

(8) 

/= E ^UA a ) (5.1) 

a=l 

and 

(ft) 

^EWA) (5-2) 

a=l 

where the positive coefficients A a are given by: 

\ a = Tr(i)A) (5.3) 
and the functions fj'(B) and gy(-B) are defined as: 

fAB) = h-r-Tr(J?,BJ?& + j*,Bjy& + JlBj*&) 

3 Ti(BW) 3 3 3 3 3 3 ' 

g r (B) = ^^Tr( (5.4) 

Here Jp, J!-,, and J?, denote the angular momentum operators truncated to the subspace 
TCjt. As before, the values of (fj a (A a ), gj a (A a )) are insensitive to the rescaling of A a by a 
multiplicative factor, so consequently X a are free variables up to the constraint: 



12 



£ X a = 1. (5.5) 

a=l 

Hence again each (/,<?) point is a convex combination of independent (/j Q (A») ) fl'ja(^a))- 
Consequently, it is sufficient to examine the much simpler case when all the X a are all 
zero except for one value of a. This means that out problem reduces to finding the upper 
boundaries of individual (fj>,gj>) regions, where f G j, j — 1, . . . , j — [j\ . 

Henceforth we drop the index a and consider a single operator Ay mapping the fully 
symmetric subspace TCj onto a certain subspace TCj> with the total angular momentum equal 
to j'. We can also assume with no loss of generality that the operator Aj> is normalized in 
such a way that: 

1 = Tr(A],A r ) 

= E E \(jW\A f \j;m)\ 2 . (5.6) 

m'=—j' m=—j 

The explicit expressions for the functions fj>(Af) and gji(Ajt) take then the following form: 

j' j 

ff( A j') = E E rn'm\(f;m'\A f \j;m)\ 2 

m'=-—j' m=—j 

j'-l 3-1 

+ E E \J(j' + m' + l)(f -m')(j +m + l)(j -m) 

m'=—j' m=—j 

xRe[(j"; m'\A f \j; m)((j"; w! + l\A f \j; m + 1})*} (5.7) 

and 

3' j 

9f(Af)= E E m\(f;m'\A f \j;m)\ 2 . (5.8) 

m'=—j' m=—j 

As the next step to simplify the problem we note that the phases of the matrix elements 
(j';m'\Aji\j;m) can be set to make all of them real and nonnegative. This maximizes the 
second term in Eq. ( |5.7| ) while leaving unchanged the expressions for gf(Aj') and Tr(A-,A,v). 

We will now demonstrate that among all the curves bounding the allowed regions of 
(fj',gj>), the curve for j' = j encompasses the largest region in the (f,g) plane, which 
includes all other bounds obtained for j' < j. For this purpose, we will show that given 
an arbitrary operator Aj> : TCj — > Tij' satisfying the condition gj'(Aji) > 0, it is possible to 
construct an operator A' : TCj — > TCj mapping the fully symmetric space such that: 

> f f (Aj,) (5.9a) 
g {A') = g r {A r ) (5.9b) 

Tr(A' f A') = Tr{A),Aj,) (5.9c) 
Hence the operator A! will be always more optimal that the original operator Aji. 
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The explicit construction of the operator A' is given by: 
(j;n\A'\j;m) = 



if; n-j + j'\A f \j; m), if n > j - 2f 
0, if n < j - 2j' 



(5.10) 



It is straightforward to check that the operator A' defined above automatically satisfies condi- 
tions given by Eqs. (|5.9b|) and ( |5.9c|) . In order to prove that condition ( p.9a|) is also satisfied, 
let us express fj{A') in terms of the matrix elements of the operator (j'; m'\Ajt\j; m): 



fM')= E E K +j- /HO*'; ™'\A f \j; m )\' 

m'=—j' m=—j 

i'-i 3-1 



+ E E \/ (2j - f + m' + l)(f - m')(j + m + - m 

m'=—j' m=—j 

xRe[(j"; m'\A r \j; m){(f] m + l\A r \j;m + 1))*} 



The second term of the above formula majorizes the second term of Eq. 
we have 



2j - f + m' + 1 > Jf + m' + 1 



(5.11) 
since for j > j' 

(5.12) 



and all the other factors have been assumed to be nonnegative. This observation can be 
combined with the decomposition of the first term in Eq. ( |5.11 ) into two parts proportional 
to m' and j — j', which yields: 



(5.13) 



Since j > j' and we have assumed that g(Ay) > 0, this proves that Eq. ( |5.9a| ) is indeed 
satisfied. Of course, the condition g(Aj>) > is fulfilled automatically for all the operators 
relevant to the trade-off, as according to our discussion from Sec. [IV] the whole region of 
interest for g is < g < j. 



VI. OPTIMIZATION 



We will now show that the search for the trade-off curve can be decomposed into a set 
of even simpler independent constrained optimization problems. To proceed further, it will 
be convenient to introduce vector notation. Let us define: 

Ik = —j + max(0, k) 

Uk = j + min(A;, 0) (6.1) 

where the index k is from the range — 2j < k < 2j. For brevity, we also denote: 

a k rn = {y,m-k\A j \y,m) (6.2) 

where Ik < m < u^., and we assume that all the matrix elements are real and nonnegative. 
We can now introduce 4j + 1 real vectors: 
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a* = (<,a? fc+1 ,...,a* fc ). (6.3) 

The length of the vector with the index k is equal to 2j + 1 — \k\. These vectors are diagonal 
stripes of the matrix (j; m\Aj\j; n). Using the vector notation, we have: 

/= E /V) (6.4a) 

k=-2j 

g= J2 */V) ( 6 - 4b ) 

k=-2j 

1= J2 h k (a k ) (6.4c) 

k=-2j 



with 



and 



u k u k — 1 

/V) = E ™(™- *)(<&) 2 + E tMA-i ( 6 - 5a ) 

ra=Z fe m=l k 

g k (a k ) = E m(a k J 2 (6.5b) 

^(a fc ) = E (6.5c) 



7, 



^/ (j — + m + l)(j + fc — m)(j — fc + m + 1). (6.6) 



We can now use the same reasoning as before to restrict our interest to a single vector a fc 
with a fixed value of k. This vector should be normalized to unity, (a k ) T a k = 1. The allowed 
region for / and g is defined as a union of regions bounded by curves obtained for different 
fc's, with the index k running from — 2j to 2j. We will now discuss several properties of 
the curves depending on the sign of k, which will allow us to restrict our search for the 
optimality curve to a smaller set of fc's. 



A. Case k < 

For negative k, we note that given an arbitrary vector a k , one can use its elements in 
the same order to construct a vector ar k simply by taking ar k = a k . A simple calculation 
shows that: 

f~ k (a' k ) = f k (si k ) 

g- k ( a ~ k ) = «?V) - k (6.7) 

Thus, one can obtain from any trade-off curve for k < a certain trade-off curve for — k > 
which is shifted along the g axis towards higher values. All the trade-off curves for k < 
are hence suboptimal, and we can further restrict our attention only to the case of k > 0. 
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B. Case k > 

We will now show that the trade-off curves obtained for k's greater or equal to y/2j lie 
completely within the region bounded by the curve corresponding to k = 0. This will allow 
us to exclude all k > ^/2j from further analysis. In order to prove the above lemma, we will 
demonstrate that the maximum value of / attained by the trade-off curves for k > y/2j lies 
below the minimum value of / on the trade-off curve obtained for k — 0. 

We start from the observation that the complete trade-off curve for k = lies above the 
value / = j 2 . Indeed, let us define the vector a = (sin%, 0, . . . , 0, cosx) with x from the 
range to 7r/4. It is straightforward to check that we have /°(a°) > j 2 over this range of 
X, whereas g°(sP) = jcos2x, which can assume any value between and j. Thus for any 
g from the range < g < j relevant to the trade-off curve we have a vector such that the 
corresponding value of / is larger or equal to j 2 . Consequently, the complete trade-off curve 
occupies the region of the (/, g) plane defined by the condition / > j 2 . 

Next, we prove in Appendix y that for an arbitrary normalized vector a fc , the function 
f k (sL k ) is bounded by: 

/V)<J'(j + l)-y- (6-8) 

The right-hand side of the above bound can be compared with the the minimum value of 
f° — j 2 on the trade-off curve for k = 0. If for a given k the general upper bound on f k given 
by the right-hand side of Eq. (|6.8|) is below the value j 2 , then the trade-off curve obtained 
for this specific k will definitely be majorized by the trade-off curve corresponding to k = 
over the whole region of interest. Hence we can exclude all k's satisfying j(j + l) — k 2 /2 < j 2 , 
which after simplification yields k > y/2j. Consequently, we can restrict our attention to 
nonnegative k's from the range: 

< k < sj2j. (6.9) 

Recalling that 2j = N, it is thus sufficient to consider k's from a finite set of only [VW] 
values. The number of independent real variables that have to be optimized for a given k is 
equal to N + 1 — k. 



VII. NUMERICAL PROCEDURE 

Our task is now reduced to finding the trade-off curves for a set of k's defined in Eq. ( |6.9| ). 
To complete this task, we shall resort to numerical means. For a given k, define {uu — h + 
1) x (v,k — h + 1) real symmetric matrices JF fc and Q k \ the matrix T k has lk(h — k), (Ik + 
l)(lk — k+ 1), . . . , Uk{uk — k) on the diagonal, and 7^/2, li k+ i, ■ ■ ■ , lu r -i on either side of the 
diagonal. The matrix Q k has Ik, h + 1, • • • , Uk on the diagonal. We can now use the method 
of Lagrange multipliers to find the maximum of / having fixed the value of g. Specifically, 
we need to optimize the expression 

f k {a k ) + \g k {si k ) - (ih k {sL k ) = (a fc ) T (^ fc + \Q k - /il)a fe (7.1) 

with the constraints: 
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g\* k ) = g 
h k (a k ) = 1 



(7.2) 



and A,/i being the Lagrange multipliers. 

Differentiating the right hand side of Eq. ( [7.1|) with respect to the elements of the vector 
a fc we obtain that the maximum occurs when 



{T k + XG k - fxl)a k = 



(7.3) 



that is, a k is an eigenvector of the matrix T k + \Q k corresponding to its maximum eigenvalue. 
Assuming that this eigenvector is normalized to one, the value of g is given by the product 
(a k ) T Q k a k , which implicitly depends on A through the vector a k . In order to plot the 
trade-off curve as a function f k (g k ) we would need to invert this relation. However, we 
can equivalently consider the trade-off curve as parameterized by the Lagrange multiplier A 
running from to oo. The value A = corresponds to optimizing f k without any constraint 
on g k , whereas in the limit A — > oo the optimization gives the maximum attainable value 
of g k and the corresponding value of f k . Thus the procedure of calculating the trade-off 
curve for a given k can be summarized as follows: for A's from the range < A < oo find 
the normalized eigenvector of T k + \Q k corresponding to the largest eigenvalue, and use it 
to calculate [a k ) T J rk a k and (a k ) T Q k a k which give the point (/ fc (A), g k (\)) on the trade-off 



We have performed numerically the task of finding the eigenvectors a k . In numerical 
calculations, it is convenient to change the parameterization of the trade-off curve according 
to A = x/(l — x) with x running from to 1, and to diagonalize the rescaled matrix 
(1 — x)T + xQ k . We note that this matrix is of the tridiagonal form, for which there exist 
efficient numerical algorithms [[HJ. After finding the normalized eigenvectors we can plot 
the parameterized trade-off curve (f k (x), g k (x)) for all relevant fc's from the range given by 
Eq. ( |6.9|) , and next find on this basis the region of fidelities allowed in quantum mechanics. 
To illustrate this procedure, we plot in Fig. |2| the trade-off curves for N = 10 and the relevant 
range of fc's. It is seen that the quantum mechanically allowed region is bounded by a single 
trade-off curve obtained for k — 0. This was the case also in all other cases we investigated 
numerically. The actual bounds on the operation fidelity F versus the estimation fidelity G 
for several values of N are shown in Fig. [|. 

On the basis of our numerical studies we conjecture that the trade-off curve for k = 
is always optimal. This conjecture is supported by the fact that both the extreme points 
corresponding to x = and x = 1 are attained only for k — 0. Indeed, we have seen in 
Sec. |VI B| that for a given nonnegative k the whole trade-off curve lies below j(j + 1) — k 2 /2. 
Consequently, the maximum value of / = j(j + 1) corresponding to unit operation fidelity 
is achieved only for k = 0. This is also the case of the other extreme point corresponding 
to the optimal estimation: it is straightforward to note that for any k > the expression 
for g k has the same maximum value equal to g = j. This value is obtained for the unique 
vector of the form a k = (0, 0, . . . , 1). The corresponding value of f k is f k = j(j — k). Thus, 
for optimized g the largest attainable value of /, equal to j 2 , is obtained only for k = 0. 



curve. 
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VIII. ATTAINABILITY OF THE BOUND 



We will now show that the trade-off curves computed in the previous section are tight, 
i.e. that they can be attained by physically realizable operations. So far, we have considered 
only the trade-off curve generated by a single operator A : 7ij — ► 'H® N satisfying the 
condition Tr^M) = 1. The critical step which allowed us to focus on a single operator was 
the replacement of the full trace-preserving condition in Eq. (fOjj) by its trace. We will now 
present a method for constructing a complete quantum operation from a single operator A 
such that it generates the same point on the fidelity trade-off curve. 

The classical outcome of the operation we construct has the continuous form of a triplet 
of Euler angles which we will denote by S. The operation element corresponding to a specific 
outcome S is given by: 



A 3 = J2j + 1U(Z)AU*(E) (8.1) 

It is straightforward to verify that the operation fidelity F and the estimation fidelity G for 
this operation are given respectively by: 

F=- + — — ^ -f(A) (8.2) 

2 2j(j + l) M } 1 ' 



and 



G = - + , 1 r g(A) (8.3) 
2 2(j + l) yv ; v ; 



where f(A) and g(A) are defined in Eqs. ( |4.8a| ) and (|4.8b ). This confirms that we obtain 



the same point on the trade-off curve as for the operator A itself. The only condition we 
need to check is the completeness of the operation on the fully symmetric subspace Jif 



I 



dEAlA s = ij. (8.4) 



Of course, outside Tij the value of this integral vanishes, as the operator A is assumed to be 
zero there and the rotation matrices i7(H) do not mix subspaces with different values of the 
angular momentum. 

In order to prove that the completeness condition given by Eq. ( |3.4|) is indeed satisfied, let 
us consider the matrix element of the left-hand side of the above expression in the eigenbasis 
of the angular momentum operator J?: 

= (2j + 1) E / d ^(r, m\U(E)\j; m')(j; m'\A^A\j; n')(j; n'\U\Z)\j; n). (8.5) 



The integral over the product of the elements of rotation matrices (J; m\U(E,)\j; m!) 
V 3 mm ,(E) can be performed explicitly using the standard formula: 
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JdZVl m ,(E)[vUm*= 5 i^- (8-6) 
With the help of the above identity we have: 
J dE (j;m\AiA s \j;n) 

= (2j + l) E JdEV^mVU^T^rn'l^A^n') 

m' ,n'=—j 

j 

= E 5 mn 5 m i n i(j]m!\A ] A\j]n) = S mn Tr(A^A) = 5 mn (8.7) 

m' ,n'=— j 

This completes the proof that the operation {A^} satisfies the full trace-preserving condition 
in the symmetric subspace of the iV-qubit Hilbert space. Consequently, the trade-off curve 
calculated for a single operator A is attained also by complete quantum operations. 



IX. DISCUSSION 

In this paper we calculated the fidelity trade-off for finite ensembles of identically pre- 
pared qubits. This trade-off relates the quality of estimating the quantum state of the qubits 
to the minimum disturbance of the original that has to be introduced in course of this pro- 
cedure. The obtained trade-off curves can also be viewed as a characterization of a specific 
asymmetric quantum cloning procedure, which given N qubits produces the same number 
of clones with a decreased fidelity F, and additionally an arbitrarily large number of clones 
with a lower fidelity G. 

The calculation of the trade-off curve was based on a combination of analytical tech- 
niques and numerical calculations. The results obtained analytically allowed us to reduce 
significantly the complexity of the optimization problem. One should note that a single 
operator A : TCj — > 'H® N mapping the fully symmetric subspace onto the whole Hilbert 
space of iV qubits contains 2 N+1 (N + 1) independent real variables. If one wanted to find 
numerically the trade-off curve assuming such general form of the operator A, the number 
of parameters in the optimization problem would explode exponentially with the size of the 
ensemble. Fortunately we were able to demonstrate that the problem of finding the trade-off 
curve can be reduced to [viV] independent optimization problems, each involving only no 
more than N + 1 real parameters. This is a substantial reduction of the problem, which 
enables one to handle numerically much larger ensembles of qubits. 

There are several elements of our work which could be investigated further. First, it 
would be interesting to prove our conjecture that the trade-off curve obtained for k = 
is always optimal. This would reduce further the complexity of the problem remaining to 
solve numerically. We have made several observations that might be helpful in this proof. 
First, numerical calculations suggest that the eigenvalues of the matrix T k considered in 
Sec. |VII| belong to the analytically defined set {—v{y — \)/2 J r 2ju—j 2 \u = 0, 1, — 
The largest of these eigenvalues is j(j + 1) — k{k + l)/2, which itself improves the upper 
bound given in Eq. (|6.8| ). Inspection of numerical results suggests also that the value of 



g (a ) corresponding to maximized / (a ) is equal to k/2. Thus, both the extreme points 
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of all trade-off curves for k > lie beneath the one obtained for k = 0. This observation, 
combined with a demonstration that the curves have appropriate monotonicity and convexity 
properties, might prove the universally optimal character of the k = curve. 

Another interesting direction is investigating in more detail quantum operations that 
saturate the trade-off inequality We have shown that given a single operator which generates 
the values of F and G lying on the trade-off curve, one can construct a complete quantum 
operation that satisfies the full trace-preserving condition. The described operation had a 
continuous classical outcome in the form of a triplet of Euler angles. It would be interesting 
to investigate operations with a finite [0] (and possibly minimal ||) number of outcomes 
that also saturate the quantum mechanical bound on the fidelities. 
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APPENDIX A: MULTIPLICITIES OF ANGULAR MOMENTUM 

REPRESENTATIONS 

In this Appendix, we derive the mutliplicities fij> of subspaces with the fixed value of the 
angular momentum j' appearing in the decomposition of the Hilbert space of N qubits. Let 
us consider the operator: 

N 

Z{(3) = (g)exp(^) = exp(20J*), (Al) 
i=i 

where the subscript i enumerates the qubits. From the tensor-product representation given 
on the left-hand side of the above formula we immediately obtain that 

Tr[Z(p)] = (2 cosh/?)". (A2) 

On the other hand, summation of the trace of the operator exp(2/3J z ) in all the subspaces 
yields: 

^)KE ,E«-^t ^ (A3) 

Comparing equal powers of in this expression with the expansion of the left-hand side of 
Eq. flA~2|) given by (2 cosh (3) N yields: 

I 2] \ I 2] \ 2? + 1 (2j + 1\ 

[j-ri •/•!.) :>,-i{j rr (A4) 

It is seen that the total number of subspaces in the decomposition ( [3.1|) is X)jv =j -_ui /V — 
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APPENDIX B: EVALUATION OF INTEGRALS K T 



In this Appendix we calculate explicitly the integrals K T defined in Eq. ( |3.8| ). For this 
purpose, it is convenient to switch to the angular momentum representation resulting from 
the decomposition of the tensor product Hilbert space of the N qubits into the direct sum 
of subspaces with the fixed value of the total angular momentum operator. In this repre- 
sentation, the state \Q)® N lies in the completely symmetric subspace characterized by the 
angular momentum j = N/2 and it is given by 

\n) 9 » = u{n)\j;j), (bi) 

where \j;j) is the eigenvector of J z corresponding to the eigenvalue j, and U(Q) is the 
rotation matrix for the angular momentum j. Of course, rotations rotations cannot transfer 
the state \j',j) beyond the fully symmetric subspace. Consequently, the operators K T are 
nonzero only in this subspace. In the basis of the eigenvectors of the operator J z , the matrix 
elements of these operators are given by: 

(j;m\K T \j;n) = [ cM k T (n)(r, m\U(ty\j; j\fr(n)\j; n) (B2) 



We shall use the standard notation from Ref. [Id*] to denote the matrix elements of the 
unitary rotation operators appearing in the above formula: 

(j;m\u(n)\ r , J ) = vi lJ (n) 

(j;j\u j m\r,n) = [v!^mY = (-i)"~^in-^)> (B3) 

where in the second line we have made use of the symmetry properties of the rotation matrix 
elements. 

The functions k T (Q) can also be expressed as elements of rotation matrices for the value 
of the total angular momentum equal to one: 

k^tt) = \/2£>i 10 (fi) 

fci(Q) = -V2V\ {n) (B4) 

With this notation, we can use the standard expression for the integrals of triple products 
of rotation matrix elements in terms of the Wigner 3-j symbols: 

(j;m|A ±1 |j; n ) = T V2(-1)-J J dQVi^V^^V^Q) 



and 



(r,m\k \ r ,n) = (-1)^/ dn2^(n)x>i B _ i (n)PJ J (n) 



■1)-M 3 3 nil' J nV (B6) 
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Inserting the explicit form of the 3-j symbols yields: 



and 



(j;m\K ±1 \j;n) = 5 m±1 , 



(j;m\K \j;n) = 5 r , 



(j ± n)(j T n + 1) 
(j + 1)(2j + 1) 



m 

(j + l)(2j + l) 



(B7) 



(B8) 



The expressions appearing on the right hand sides of the above formulas are proportional 
to the operators J T and J z restricted to the completely symmetric subspace with the total 
angular momentum j. Thus, the integrals K T can be conveniently written as 



and 



K ±l 



K 



(j + l)(2j 



1 



(j + l)(2j + l 



(B9) 



(BIO) 



where by the subscript j we have denoted the angular momentum operators truncated to 
the subspace characterized by the angular momentum j. 



APPENDIX C: DERIVATION OF THE BOUND ON f k (a k ^ 



Here derive an upper bound on the functions f k (a k ), defined in Eq. ( p.5a| ), for k > 0. 
For this purpose, let us rewrite 7^ defined in Eq. (|6.tj|) to the form: 



1 rn 



\ 



1 kv (e 
J~ m + 2 - 5. 



j + m + l- - 



Using the inequality: 



\J (x\ — x\) [x\ — xl) < X1X3 — X2X4 

valid for any real x\ > X2 > and X3 > 24 > 0, we obtain that 



>y k < 5 

lm — "mi 



where: 



5 k m <j(j + l)-[rn 



m hi 

2 



(Cl) 



(C2) 



(C3) 



(C4) 



We will also use the inequality between the geometric and arithmetic means to estimate: 

1 



°m a m+l — 2 t( a m) ( a m+l) ]■ 



(C5) 
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We can now use the inequalities given in Eqs. (|C3| ) and (|C5|) to find an upper bound on the 



second term in Eq. (|6.5a ). A simple rearrangement of the terms yields: 



f k (a k ) < £ m{m - k){a k m f + - £ 5 k m [{a k m f + ,„ 



k N2n 
m+l) J 



«fc i i «fc-l i 



£ m(m - ^)(a^ t ) 2 + -<(afj 2 + 2<-i«) 2 + £ + 6 ™ ){a ' ; 



(C6) 



We can now add to the right-hand side two terms of the form 5 k k _ 1 (af k ) 2 /2 and 5* (a k J" 2 /2 
which are nonnegative as can be easily checked, and use the identity: 



+ <t) = J(J + 1) - m(m - k) - y. (C7) 



This finally yields: 



u k u k / U2 N 

/V) < £ m(m - fc)(<4) 2 + £ j(j + 1) - m(m - *) - - ] (a 

m=lk iri=lk ^ ' 



fc \2 



+ X>™) 2 - (C8) 

m=lk 



If the vector a fc is normalized to unity, i.e. ElLiJii) 2 — 1> we obtain the bound given in 
Eq. (H). 

We note that for /c = the derived bound is tight as it can be achieved for a fc_0 = 
(W2TTT,...,W2J+T). 
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FIGURES 
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FIG. 1. The trade-off curves between the operation fidelity F and the operation fidelity G, 
depicted using absolute (top) and relative (bottom) scaling of the axes, for several values of the 
ensemble size: N = 1 (solid), N = 5 (dashed), N = 20 (dotted-dashed) , and N = 100 (dotted). 
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FIG. 2. The trade-off curves for the ensemble of N = 10 qubits obtained by solving the opti- 
mization problem for k = (solid), k = 1 (dashed), k = 2 (dotted-dashed), and A: = 3 (dotted). 
The allowed region is a union of the regions bounded by each of these curves. 
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